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ABSTRACT 



■ When electron-electron correlations are important, it is often necessary to use "exact" numerical 

■ methods, such as Lanczos diagonalization, to study the full many-body Hamiltonian. Unfortunately, 



such exact diagonalization methods are restricted to small system sizes. We show that if the Hubbard 
U term is replaced by a "periodic Hubbard" term, the full many body Hamiltonian may be exactly 
solved, even for very large systems, though for low fillings. However, for half-filled systems and large 
U this approach is not only no longer exact, it no longer improves extrapolation to larger systems. 
pH ■ We discuss how generalized "randomized variable averaging" (RVA) or "phase randomization" schemes 

O ■ can be reliably employed to improve extrapolation to large system sizes in this regime. This general 

approach can be combined with any many-body method and is thus of broad interest and applicability. 

X 

H ■ INTRODUCTION 



Solutions of the Peierls-Hubbard Hamiltonian (PHH) [1] for small clusters are strongly dependent 
on the boundary condition (BC); i.e., for the case of polyacetylene, whether the carbon atoms are 
viewed as being on a real ring, such as benzene, a real chain, such as (l,3,5)-hexatriene, or a more 
exotic geometry, such as anti-periodic BCs. We have found that [2,3] to reduce finite size corrections 
to the calculated optical gaps and spectra, thus improving the extrapolation to the infinite case, it is 
effective to employ novel "averaging" techniques to "randomize" the many-body levels which influence 
those correlation functions of interest. Briefly, by randomization, we mean simply that the total energy 
of the system is to be viewed as a weighted average of the energies derived separately for each of several 
different "random" values (RV) of a particular parameter, which may be the BC, a hopping integral, 
an on-site energy, or the Hubbard U: E = X]{_Ry} ^RvE[RV], where xrv is a normalized weighting 
factor with X]{_Ry} ^RV = 1- Similarly, any correlation function can be viewed as the average of of 
its values from the separate RVs. The number of RV, N^y, may be varied, and the choice of the RV 
and xrv rnay be truly random, or may be a pre-specified set of values, which may in turn depend on 
the quantity being studied. The case of phase randomization and its relation to Bloch's theorem are 
described below. More generally, this randomized-variable-averaging (RVA) technique may be viewed 
as a procedure to mimic disorder. 



We focus in this manuscript on the calculation via exact Lanczos diagonalization of the self- 
consistent uniform dimerization and one-photon optical absorption of the ground state within the 1-D, 
1/2-filled, one-band PHH. However, RVA is equally applicable to the calculation of, e.g., the phonon 
modes of the ground state, as well as the geometry and self-consistent absorptions of doped or higher 
lying {e.g., triplet) states [3], luminescence spectra, structure factors, etc., in single- or multi-band 
models in 1-, 2-, and 3-D, via Lanczos diagonalization, as well as other numerical procedures on small 
lattices, such as Monte Carlo. 

PHASE BOUNDARY CONDITION AVERAGING AND BLOCK'S THEOREM 

Bloch's theorem tells us that the single particle wavefunctions of a system of size M-N with 
periodicity N and periodic boundary conditions (PBCs) are exactly the single particle wavefunctions 

of the ensemble of systems of size N for each of the phase BCs -^{N + 1) = e'*^{l), cf) = 
t = 1,...,M. Restated, the (exact) properties of a larger (single particle) system may be found by 
forming a (particular) average over smaller systems with different BCs. This "boundary condition 
averaging" (BCA) or "phase randomization" is a special case of our general RVA procedure. Blochs 
theorem may be generalized to many particle wavefunctions and is applicable here if the form of the 
Hamiltonian is modified slightly. The one-dimensional PHH is 

If we replace the Hubbard term, U^^^ n^^^Ui^i, by its periodic analog, 

jj M N 

Hei-ei '^ne+^.N,'fne+u.N,i , (2) 

ij,,i/=i e=i 

then, if the original problem had periodicity N, {ne^cr) = {'ne+vN,a), we have, in a mean-field sense, 
the "same" Hamiltonian. With this new "periodic-Hubbard" Hamiltonian (periodic PHH) the Bloch 
analysis of a given large (periodic) system {N -M sites, N^-M electrons) is accomplished by considering 
a small system {N sites, Ne electrons) with several different BCs. The method scales linearly with M, 
allowing one to handle reasonably large systems, and sets a sound theoretical basis for the empirical 
observation [2] that RVA can be used to smooth optical absorption spectra obtained via exact finite-size 
diagonalization. 

Using the periodic-Hubbard term, Eq. (2), assuming the lattice distortion has the same periodicity, 
S{£ + N) = S{£), and using Y.^=i Y,i=i Z]o-=r,i '^^+m-JV,o- = Ne-M, one can show: 



where 



^ M N 



M,i^=l i,j,k,l=l Si,S2=T>J. 



H{i, j, k, I; si , S2) =H{i + N,j + N, k, I; 81,82) 

=H{i,j,k,l;s2,8i) = H{k,l,i,j; 81,82) = H{j,i,k,l;8i,82) 
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Note that H{i,j,k,l;si,S2) is independent of /x, u, and M. The Bloch analysis on the many-body 
eigenfunctions using the symmetries of H{i,j, k, I; 81,82) leds to a ^ of the form: 

N 

*fei,...,fcp = Yl V'fei,...,fcp(ni,-,np) 

c7i,...,<Tp=t,i n,i,...,np=l 

(4) 

Atl = l Mp = l 

where V'fci,...,/c,,...,fcp(ni, .-.jnp) = tpki,...,ke,...,kp{ni, ...,ne+fie-N, ...,np) , ke = (27rj/M), je{0, ...,M- 

1}, and p = Ne-M.We find for i,je{l, ...,Ar}, 

M N N 

J2 ^l+'^-N,sCj+C^+uo)-N,s'^ku-,kp = Y '^'^^^-^k^'^.-^np) dn^,jd„^,s e'^'''° 
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V pifciMi V ... V \ " J pikefJ-e v . . . v \^ r'^ pikpfip 

Mi=l «=1 Mp=l 

(5) 

(fo comes into play for operators like cjy ^cn+i,s)- From Eq. (5) we see Hp^r is diagonal in the fc's, 

Hper K,...,kp = , (6) 

as are p = J^^^^ ni^a- and J = i .^(-to + Oih){c\^„Ci+i,a - c]+i,crC^,cr), and the eigenfunctions of iJper 
can be written in the form Eq. (4). 

The symmetry of ^ implies that if a^^ = at^ and k^^ = ke^ then n^^ 7^ ne^ . Thus there are at most 
2A'' of the ke the same. We postulate that the ground state lies in the manifold with each of the M 
distinct values for ki occurring times and (for the half-filled band) with equal numbers of up and 
down spins (in general, the largest manifold). This can be checked, and we stress that the analysis up 
to this point is exact for the periodic PHH. In this manifold we can write ^' as: 

2N 
ai,...,aM=l 

where qi = (27r£/M) and 

N 

K;^ = E E E <ainU:;nNj 
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where A^^ = [[^^^^]], = [[^]] ([[ ]] denotes integer value of). Note 

4+^,,CA,,. ^ll, = e"^ cl^c^,^ *5 , (8) 



and so solving the electron problem on the full N-M sites with PBCs for fixed q is exactly equivalent 
to solving 

iiper ^q,a q,a ^ q,a \^ ) 

on sites with the g-dependcnt BC defined by Eq. (8). Thus, to find the exact electron eigenfunc- 
tions ^f^^ on iV • M sites, we only need to solve Eq. (9) with N sites and A^e electrons for each of the 
M BCs, rather than the full A'' • M site problem. 



To obtain the exact eigenfunctions of the full Nf,-M electron problem, we need to solve for 
^^{ai, ckm)- This is where an approximation must be made to be able to solve the problem numer- 
ically, as we assume N is already as large as computationally feasible. The full eigenvalue problem 
yields 
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(10) 



To solve Eq. (6) we must use some approximation for Eq. (10), such as perturbation theory or mean- 
field. Thus BCA means we treat electron-electron correlations within the electron manifold exactly, 
and electron-electron correlations between A^e electron manifolds approximately. The BCA results 
reported here are zero-order perturbation theory: we have assumed (j)^{ai, ...,q;m) to be a product of 
5-functions; i.e., for the ground state we use: 
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(11) 



Note that if we had been interested in systems at low filling (Ne electrons), rather than near half filling 
(A^e ■ M electrons), there would be no need to make any approximations. 



Comparison of the periodic and standard Peierls-Hubbard Models at half-filling 

At U=0 Bloch's theorem for the many particle wavefunctions as formulated above is exact at all 
fillings. For the half-filled PHH, we have just shown that phase BCA involves two approximations: 
first, replacing the Hubbard term by Eq. (2), and, second, ignoring correlations between the particle 
wavefunction leading to Eq. (11). We now test these approximations. In Fig. 1 we show the minimum 
energy dimerization of an 8-site system as a function of U for (i) phase averaged solution and (ii) 
the exact solution of the half-filled band with the periodic Hubbard term and periodicity 2,4, and 
8 (periodicity 8 is the the usual Hubbard). The approximate periodic (BCA) solution lies between 
the exact periodic solution and the exact pure-Hubbard solution. Since the original problem was for 
U ■ M J2f=]^ '^^n^i on N ■ M sites, one might argue that the approximate phase-averaged result for 
U ^iVf-ei on N sites may be closer to the desired answer than the exact periodic-Hubbard N ■ M 
site result. 

In Fig. 2 we show the dimerization amplitude obtained by phase BCA. We see that the extrapolated 
infinite behavior at small to intermediate U is well approximated after phase averaging for A''=10 even 
without extrapolation. The agreement is best for small U, while no change from the pure periodic 
behavior is found for large U. That this particular technique is expected to have no eff'ect at large U 
can be seen by examining the effective spin-Peierls Hamiltonian [4], which is asymptotically independent 
of the phase of the BC. However, the phase-averaged analysis does lead to the interpretation of the 
finite-size result as the result for the larger system with an effective periodic-Hubbard interaction. 
Thus, if U is scaled by the system size, as is done in Fig. 2, the infinite behavior should be more easily 
extrapolated. It appears that this is at least approximately true. 

To show how phase BCA affects the optical spectrum, in Fig. 3a we show the (Lorentzian broad- 
ened) spectrum at C//to=0.4 for A^=8 and 5 phase BCs (enough to yield a smooth spectrum at U=0). 
It is clear how the different BCs "fill in" the spectrum. Fig. 3b shows the spectrum at U/tQ=l.6. Here, 
despite the several BCs used, the spectrum remains sparse, due to the loss of effectiveness of the phase 
BCA scheme as discussed above. If we were to use a broad enough Lorentzian to smooth this spectrum 
(the same width as in Fig. 3a was used), the optical gap edge would be lost. 



Fig. 1 (left). Comparison of the approximate and exact periodic-Hubbard minimum energy dimerization 
as a function of U for to=2.5 eV, a=4.1 eV/A, and K=21 eV/A^. 

Fig. 2 (right). Dimerization vs. U for various N: phase BCA. Parameters as Fig. 1. 



Fig. 3. The Lorcntzian broadened optical absorption spectra calculated using the phase BCA method 
for (a) U=l eV and (b) [/=4 eV. Note the "sparseness" of the spectrum in the latter case indicating 
the failure of the phase BCA approach. For comparison, we also show the spectrum at U=10 eV using 
(c) amplitude BCA and (d) an RVA procedure where U is varied on one site. Note these RVA spectra 
are dense and show the same features. Parameters as Fig. 1 but 6=0.14 A. 



AMPLITUDE BOUNDARY CONDITION AVERAGING 

The modification of the phase of the BC discussed above can also be viewed as passing a random 
magnetic flux through an ensemble of closed ring and studying average properties. In this sense we have 
"randomized" the locations of the momentum space states. We have already indicated that RVA can 
help us transcend the limitations of this "bond phase/magnetic flux" approach. We can "randomize" 
electronic properties by, e.g., changing a local hopping or an on-site energy or a Coulomb repulsion 
somewhere on the chain. One could also introduce an additional field and vary it about zero. Although 
there has been some limited earlier work on using modified BCs in the context of the Hubbard model 



[5], there is as yet no provably accurate prescription for for arbitrary U and V. However, certain 
intuitive rules must guide us. First, whatever change is made to randomize must, of course, do so 
effectively. Second, the change in the system must be negligible as the lattice size is increased to 
infinity; for example, if only one bond or site is varied from calculation to calculation, then the effect 
of such a change is immaterial in the thermodynamic limit. Finally, the behavior for small lattice 
sizes must be illustrative of the infinite-size limit. Put another way, one must still be able to make a 
reasonable extrapolation to the infinite chain. 

We have found [2] an "amplitude BCA" ("scaled-hopping") technique to be effective and pro- 
duce results in good agreement with expectations based on both strong- and weak-coupling analytic 
arguments. For this method, we randomize by varying the magnitude rather than the phase of the 
"boundary" hopping - i.e., between sites 1 and A'' - typically from —to to +to in ten to twenty equal 
steps. The individual spectra are then added together with weights Xi chosen to minimize the total 
length of the final curve, though giving each spectra equal weight yields virtually identical results. This 
clearly incorporates the special case of using only periodic {x=l) or antipcriodic (x= — 1) rings, the 
JT/nJT difference being important for weak coupling. It also incorporates the case of the open chain 
{x=0), which has the "best" single BC size dependence (though with attendant "end effect" problems). 
Finally, it works in the strong-coupling limit. 

In Fig. 3c, we show the spectra produced by the amplitude BCA technique for for a larger U than 
that in Fig. 3b where phase BCA failed. In Ref. [2] we showed that not only are the gross features 
of the spectra obtained by amplitude BCA in agreement with strong-coupling calculations [6], but in 
addition they show substantially more interesting detail, such as the "decoupled-dimer" peak located 
at ~ U/2+ \J {U /2)'^ + (2to)^ in systems with strong electron-phonon and electron-electron couplings. 
To emphasize the generality, in Fig. 3d we show the spectrum obtained by varying U on the first site 
between and 20 eV. This large U variation produces states in the optical gap, but otherwise the 
spectrum is unchanged from Fig. 3c. We feel that these results, and those in Ref.s [2,3], demonstrate 
that RVA is an effective means for reducing (and in some cases practically eliminating) finite size 
dependence, yielding results that can be confidently extrapolated to the infinite size limit, and which 
are in good agreement with known analytic results. 
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